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ABSTRACT 

We present a study of the submillimeter (submm) emission of Distant Red Galaxies (DRGs) . The 
DRGs are selected by the criterion J — K > 2.3, and are generally massive galaxies at redshifts higher 
than 2, with red rest-frame optical colors. Using a deep SCUBA submm image of a field centred on 
the cluster MS 1054—03, we obtain a statistical detection of the DRGs at redshift z = 2 — 3.5, with 
an average 850 /im flux density of 1.11 ± 0.28 mJy. The detection implies an average star formation 
rate (SFR) of 127 ± 34M0yr -1 (lensing corrected), assuming that the far-infrared (FIR) spectral 
energy distribution (SED) is well-described by a modified blackbody. The SFR derived from the 
submm agrees well with SFRs derived from SED fitting of optical-near-infrared data and average X- 
ray emission. Constant Star Formation models imply ages of 2Gyr, extinction Ay = 2.4 mag, which 
is consistent with the FIR to rest- frame optical luminosity ratio of ~ 15. DRGs are older and have 
lower SFRs relative to optical luminosity than (ultra-) luminous infrared galaxies, although their FIR 
luminosities are similar. The DRGs at 2 < z < 3.5 and the Extremely Red Objects (7 — K > 4) at 
1 < z < 2, which were also investigated, contribute 5.7 and 5.9 Jy deg -2 respectively to the submm 
background. Simple estimates suggest that these populations contribute ~ 50% of the flux from 
sources with 0.5 < /sso < 5 mJy, which is where the peak of energy is produced. We have therefore 
uncovered one of the most important populations of galaxies contributing to the sub-mm background. 
Subject headings: Galaxies: high- redshift — galaxies: evolution — submillimeter 



1. INTRODUCTION 

Members of the recently defined galaxy population 
Distant Red Galaxies (DRGs) are identified by their red 
observed near-infrared col ors, J — K > 2.3 (iFraruL^LaLl 
20031 Ivan Dokkum et all l200l iForster Schreiber et all 
20041 Thenceforth FS04)). Galaxies selected by this color 
cut are generally at z ~ 2 — 3.5, and are hence selected 
on their red colors in the rest-frame optical. Their stel- 
lar populations are relatively old, or are highly dust ob- 
scured. DRGs account for over half of the cosmic stel- 
lar mass at the high mass end at redshifts 2 < z < 3.5, 
somewhat more th an the contribution fr om Lyman Break 
Galaxies (LBGs) (IRudnick et al.ll20C)3t FS04). One of 
the most pressing questions concerning these galaxies is 
how high their star formation rates (SFRs) are. Based 
on near-infrared (near-IR) spectroscopy, high star forma- 
tion rates of 200 — 4OOM yr~ 1 and hig h stellar masses 
of (1 —5) x 10 11 M w were determined l)van Dokkum et all 
2004). Unfortunately, near-IR spectroscopy is available 
for only a few galaxies at this moment. Star formation 
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rates can also be estimated using other techniques, espe- 
cially using radio, X-ray or submillimeter (submm) ob- 
servations. Most of the DRGs are too faint to be de- 
tected individually at these wavelengths, but it is pos- 
sible to estimate their average star formation rate by 
"stacking" the fluxe s for a large number of sources. Us- 
ing this technique, iRubin et alJ l)2004|) have detected 
the X-ray emission of DRGs in the field of the cluster 
MS 1054—03 and derived an average star formation rate 
of 214 ± OQMeyr- 1 if all the X-ray flux is attributed 
to the SFR. However, the flux contribution from faint 
AGNs remains unclear. 

Here we present a different SFR measurement for 
the DRGs, based on the average 850 /im flux density, 
which we attribute to emission from dust that has 
been heated by young stars. We have taken very deep 
SCUBA 850 u m. observations of the field of the cluster 
MS 1054-03 l)Knudsenl 12004 Knudsen et al. in prep). 
These deep data provide a unique opportunity for study- 
ing the submm properties of DRGs. Throughout this 
Letter we will assume an ft — 0.3, A = 0.7 cosmology 
with Hq = 70kms _1 Mpc -1 . Unless otherwise stated, 
we use Vega magnitudes. 

2. SUBMILLIMETER DATA 

Submillimeter data for MS 1054—03 were obtained 
using the Submill imetre Common- User Bolometer Ar- 
ray (SCUBA; iHolland et al] Il999f) at 850 ^m on the 
15 m James Clerk Maxwell Telescope at Mauna Kea, 
Hawaii. Total integration time was 49 hours, distributed 
almost evenly over 3 pointing centers to cover the en- 
tire fiel d of interest. The data were reduced using 
SURF Ij.lenness Rr. Lightfootl ITOfll . and sources were 
extracted using the Mexican Hat Wavelet algorithm 
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(jBarnard et al.ll2004t IKnudsen et alJl2005|) . Noise prop- 
erties of the data were determined using Monte Carlo 
simulations of artificial datasets. Details of these re- 
ductio n an d analysis methods ar e described in IKnudsenl 
( 2004) and IKnudsen et alJ {2005). The total area of the 
map is 14.4arcmin 2 , though for the analysis the outer 
23" edge has been trimmed off due to high noise leaving 
a useful area is 12.8arcmin 2 . The image has an area- 
weighted la r.m.s. noise of ~ 1.2 mJy beam -1 , with vari- 
ations between 0.9 and 1.6mJybeam _1 . The flux cali- 
bration uncertainty at 850 /im is about 10%. The Full 
Width Half Maximum (FWHM) of a point source in the 
map is 15". The astrometric accuracy is approximately 
3", due to the pointing uncertainty of the JCMT. 

Nine submm sources are detected with flux densities of 
3.5 — 5.0 mJy with signal-to-noise ratio > 3. Two of the 
nine submm sources have been securely identified, both 
with DRGs, namely the galax ies MS-1383 which has a 
spectroscopic redshift of 2.423 l|van Dokkum et aLll2004l) 
and MS-723 which has a photometric redshift of 1. 88^02 
l)Forster Schreiber et al.ll2005l) . Five of the other submm 
sources have likely identifications, with at least one DRG 
or Extremely Red Object (ERO, defined by an observed 
color I — K > 4) nearby. The two remaining sources have 
no obvio us counterpart. A mas s model for the foreground 
cluster (iHoekstra et al.l l2000|) shows that the gravita- 
tional lensing magnification is small for the sources, on 
average around 20%. 

3. SUBMM EMISSION FROM DISTANT RED GALAXIES 
3.1. Statistical analysis 

We measured the characteristic submm fluxes for 
DRGs as a class of galaxies by measuring the flux in the 
850 /im map at the position of all the galaxies in the sam- 
ple. We found an average 850 /^m flux of 1.19 ± 0.22 mJy 
for 30 DRGs with K < 22.5, and with an uncertainty on 
their if -band magnitude smaller than 0.1, and minimum 
weight in the photometry > 0.2 (this measure is related 
to the total integration time - see FS04), present within 
the area covered by the SCUBA map. As the sensitivity 
is not uniform across the field, we determine the weighted 
mean and standard deviation, weighted by the squared 
reciprocal of the la r.m.s. noise at the position. The 
weighted mean compared to the unweighted mean only 
changes the result by a few percent, indicating that the 
weighted result is not dominated by a few data points 
with high weig hts. iWebb et all J2004) discussed a pos- 
sible bias in this procedure due to confusion. If two or 
more DRGs are separated by less than the 850 /im beam, 
emission might be counted twice. This is corrected for by 
estimating the flux contribution from each of the close- 
by DRGs to the flux measured at the given position (for 
details see: IWebb et al.ll2004fl . We find that this effect 
influences our measurements by ~ 2—10%. 

As mentioned above, two DRGs are identified as the 
counterparts of two individual SCUBA detections. Ad- 
ditionally, three other DRGs are found within the identi- 
fication search radius of three other SCUBA detections, 
though these are not secure identifications. We calculate 
the mean 850 /im flux density, both including and exclud- 
ing these five objects. Results are summarized in Table 

If we omit the DRGs associated with discrete SCUBA 
sources, we obtain an average flux of 0.74 ± 0.24 mJy, 
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Fig. 1. — Postage stamp sized sub-images from the sta- 
tistical stacking analysis of the DRGs in the field MS1054— 03. 
Each sub-image is 50" X 50". The images represent a weighted 
mean. The overlayed contours represent a mean 850 (im flux of 
0.0,0.2,. .,1.2mJy. Top row: DRGs for all redshifts, including (left) 
and excluding (right) sources associated with the discrete SCUBA 
sources. Bottom row: sources at redshift 2 — 3.5, left are DRGs, 
right non-DRGs. 



demonstrating that the detection is not caused solely by 
these galaxies. We show the average images of both sam- 
ples in Fig. n The emission is well-centered. 

For DRGs restricted to the redshift interval 2 < z < 3.5, 
we obtain an average flux of 1.11 ± 0.28 mJy for the 
DRGs. For comparison, the average flux for galaxies at 
2 < z < 3.5 which are not DRGs, i.e. J — K < 2.3 (hence- 
forth refered to as non-DRGs), is substantially lower at 
only 0.14 ±0.15. 

Following FS04, we fitted the full optical/near-IR spec- 
tral energy distribution (SED) with stellar population 
models for DRGs with K <21.7 and 2<z<3.5. The 
sample was separated into those best fit by constant 
star formation models (CSF) or single age stellar popula- 
tions (SSP). The CSF-subsample has an average flux of 
1.67±0.38mJy, while the SSP-subsample has an average 
flux of 0.67 ± 0.54. Excluding the DRGs associated with 
discrete SCUBA sources affects only the CSF subsample 
and its mean flux decreased to 1.15 ± 0.43 mJy. We note 
that the lensing correction is less than or comparable to 
the uncertainty of the average fluxes. 

3.2. Monte Carlo simulations 

To test the reliability of the detection of the mean 
850 /im flux from the DRGs, Monte Carlo (MC) simu- 
lations were performed by making the statistical mea- 
surements at random positions in the field. The MC 
simulations were repeated 10000 times. As with the real 
data, the MC measurements were corrected for the con- 
fusion of two or more sources separated by less than the 
distance of a beam. For 45-50 random positions on the 
map, the weighted mean is 0.13±0.17mJy. As the noise 
distribution deviates from Gaussian, the probability for 
getting a > 3u detection at random positions is 2.4%. 
The probabilities for a > 4a detection and a > 5a detec- 
tion are 0.2% and 0.02% respectively. A similar result 
is seen for 20-30 and for 10-20 random positions. Thus 
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TABLE 1 

The mean 850 /im flux of DRGs 



Including all Excluding discrete sources 



Sample 
(1) 


^ limit 
(2) 


N 
(3) 


{/sso) 
(4) 


% EBL 

(5) 


(SFR) 
(6) 


(Lv) 
(7) 


(SFR)/(L V ) 
(8) 


N 
(9) 


(/S50> 

(10) 


(SFR)/(L V ) 
(11) 


all-z 

2<z<3.5 
no-DRGs 
2<z<3.5 

CSF 
SSP 


22.5 
22.5 

22.5 

21.7 
21.7 


30 
18 

71 

9 

5 


1.19 ±0.22 
1.11 ±0.28 

0.14 ±0.15 

1.67 ±0.38 
0.67 ±0.54 


21.8(30.9) 
13.0(18.4) 

6.5 (9.3) 

10.3(14.6) 

1.6 (2.3) 


169 
159 

20 

238 
96 


7.3 
7.9 

4.7 

10.4 
6.1 


23.2 
20.2 

4.3 

22.9 
15.7 


25 
15 

67 
7 


0.74 ±0.24 
0.69 ±0.30 

0.04 ±0.15 
1.15 ±0.43 


15.1 
13.8 

1.2 
17.7 



Uncorrected for the average gravitational lensing of 20%. ^^^^^^^^ 

Units: (2) [mag]; (4) [mjy]; (5) % of the EBL value from lFixsen et alj <1998l) IPuget et"afl <1996T) 1: (6) [M Q /yr] ; 
(7) [10 10 L e ]; (8) [M /yr/lO lo L ]; (10) [mjy]; (11) [M /yr/lO 10 L ] . 



we conclude that the detected mean 850 /xm flux is real. 

3.3. Possible systematic uncertainties 

We performed simulations to verify how reliable the 
average detections were. We assumed that the star for- 
mation rate is proportional to restframe Ly luminosity, 
with different ratios for different classes of objects. We 
found that systematic errors occur at the level of 0.1-0.2 
mjy for the average fluxes. For the DRGs and EROs 
this is not a concern, but for the fainter objects (e.g., the 
non-DRGs with average flux 0.14 mjy) this implies that 
the systematic errors are at least as large as the detec- 
tions. The systematic problems are due to a combination 
of effects, including the negative "sidelobes" in the beam 
pattern, which cause the total flux of the map to be zero, 
and can reduce the flux for the faint, but abundant non- 
DRGs, and superpositions of objects. Obviously, simi- 
lar observations of other fields would help to constrain 
the errors for the DRGs; but it is unclear whether such 
measurements will ever constrain the flux from the faint 
non-DRGs. 

4. DISCUSSION 

4.1. Implied star formation rates of the DRGs 

The average SFR of the DRGs can be estimated from 
the average submm flux assuming an SED and an ini- 
tial mass function (IMF). As the submm and far-infrared 
(FIR) SED of the DRGs is unknown, we u se the SED de- 
script ion for dusty starburst galaxies from lYun fc CariTTil 
(2002) for calculating the FIR luminosity. Because of 
the large negative fc-correction in the submm, the ob- 
served 850 /im flux density is essentially constant be- 
tween redshift 1 and 8 for a given luminosity. For the 
DRG 2<z<3.5 sample, the average FIR luminosity 
is 1.2 x 10 12 L©, which is comparable to the luminos- 
ity of the local ultraluminous infrared galaxy (ULIRG) 
Arp220. Furthermore, we ba se the conv e rsion between 
FIR luminosity and SFR on iKennicuttl l)1998|) . where 
a Salpeter IMF in the mass range 0. 1-100 M© is as- 
sumed. However, we assumed a mean age of the stellar 
population of lGyr and a constant star formation dur- 
ing that time, which lowers the conversion by a factor 
of 1.5: 5F^[M yr- 1 ]= 1.17 x W- w L FIR [h & }. We as- 
sumed this longer mean age since modeling of the optical- 
Near IR SEDs suggest ages of 1-2 Gyr (FS04). We ob- 



tain an average SFR of 159±42M yr _1 for 2 < z < 3.5, 
where the error reflects the accuracy of the submm de- 
tection. When correcting for gravitational lensing, this 
is an (SFR) of 127 ± 34M yr" 1 . We note that the FIR 
luminosity and SFR have considerable uncertainties as 
they depend on the assumed shape of the SED and IMF. 

The average SFR derived from the thermal dust emis- 
sion in the far-IR (X rest ~ 250 /im) can be compared to 
the SFRs inferred from the SED fits to the optical and 
near-IR data (FS04). The SFRs from SED fits are rather 
model-dependent. One of the main uncertainties is the 
parametrization of the star formation history (constant, 
declining etc). The maximum is generally given by con- 
stant formation rates, and is on average ^170MQyr _1 
for the sample studied here (FS04). This compares well 
with the SFR derived from the submm flux. 

The good agreement suggests that the majority of 
DRGs are well approximated by simple models with con- 
stant star formation rate, ages of 2 Gyr, and high ex- 
tinction Ay — 2.4 (FS04). Despite the abundant dust 
the overall ratio of Lfir/ L opt ~ 15, which is in good 
agreement with Ay, is lower than that of ULIRGs in 
the nearby universe. This implies that DRGs are form- 
ing stars more steadily over longer periodes than nearby 
starbursts. We expect to find very large reservoirs of cold 
gas to sustain these high star formation rates. 

When interpreting the average rates it should be noted 
that not all DRGs are likely to be formin g stars. In a 
deep S pitzer study of DRGs in HDF-South, Labb e et alJ 
( 2005) found that 30% are best fit by old stellar popu- 
lations at rcdshifts higher than 2, and the remaining are 
best fit by constant star formation models. The results 
suggest that at redshifts higher than 2, the majority of 
massive galaxies were still forming stars at rather high 
rates, enough to build up their stellar mass in 1-2 Gyr. 

4.2. Contribution to the comoving star formation 
density 

We derive a comoving star formation density of about 
O^Mgyr^Mpc- 3 for the DRGs with 2 < z < 3.5. 
This can be compared with the comoving star formation 
density derived from UV selected galaxies at z = 2 — 3.5 
of approximately 0.022 to 0.11 M yr _1 Mpc" 3 (Steidel 
et al 1999), excluding, and including a correction for dust 
absorption. As the dust correction is uncertain for UV 
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selected galaxies, especially at the faint end which con- 
tributes a large fraction of the flux, we have to conclude 
that the contributions of the two populations are similar, 
but still uncertain. The main uncertainty for the DRGs 
is the luminosity function at the faint end, and the main 
uncertainty for the LBGs is the dust correction. 

We note that several authors have tried to detect the 
average submm flux of the LBGs (e.g., Chapman et al. 
2000, Webb et al 2003), without statistically significant 
results. Since our simulations suggest that systematic 
effects can be large at the very faint end, the average 
850 /im flux of these sources remains an open question. 

4.3. Dissecting the contribution to the sub-mm 
background 

Each population of high redshift galaxies contributes 
to the submm Extragalactic Background Light (EBL). 
Previous studies have addressed this for the LBG and 
the ERO populations. There has been no significant 
statistical detection of the subm m emission from LBGs. 
Through SCUBA observations, IChapman et~aTl (j2000f) 
placed upper limits on the LBG contribution to the 
submm EBL of 0.2%, while iWebb et alJ pio^ placed 
an upper limit of 20% when extrapolating to a redshift 
interval 1 < z < 5. Contrary to the LBG case, statis- 
tical detection of the ERO p opulation has been suc- 
cessful. IWehner et all l|2002f) and IWebb et alJ pOOl 
find an average submm flux of 1.58 ± 0.13mJy and 
0.40 ± 0.07 mJy respectively, and determine the contri- 
butions to the submm EBL to be about half and 10%. 
The discrepancy between the results can be primarily 
assigned to the difference in depth of the ERO selection 
an d the strong clusterin g of EROs, as the fields studied 
by Wchnc r et alJ ((2002) appear to have a large overden- 
sity of EROs. 

Using the present study we constrain the contribution 
of the DRGs to the submm EBL. For the whole sam- 
ple of DRGs in the MS 1054-03 field, the total flux 
density is approximately 9.6 Jy /deg 2 , which is 21 .8% 
of the EBL value measured bv IFixsen et al. | 119981) of 
44 Jy/ deg 2 and 30.9% of the value from iPuget et alJ 
(1996) of 31 Jy/deg 2 . For the DRG sample limited 
to the redshift range 2 < z < 3.5 the contribution 
is 5.7 Jy/deg 2 i.e. 12.9% [18.4%J of the result from 
IFixsen et alTl|1998ft jPuget et alJ l|1996|) ]. We can per- 
form the same measurement for the EROs in the interval 
1 < z < 2. We apply similar magnitude and uncertainty 
cuts as for the DRGs (see section OJ. We obtain an av- 
erage flux of 0.69 ± 0.21 mJy, by averaging 35 sources, 
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which implies a comoving star formation rate density of 
O.OSSMQyr^Mpc- 3 . The total contribution from the 
EROs is 5.9 Jy/deg 2 . Hence the total contribution to the 
submm EBL of the EROs at 1 < z < 2 and the DRGs at 
2 < z < 3.5 is 11.6 Jy/deg 2 . 

The surface density of the two populations combined 
is ~ 5arcmin -2 down to K = 22.5 mag, similar to the 
surface density derived from the submm number counts, 
where 7V(>0.5mJy) ~ 6.5arcmin~ 2 . 

We can use our measured ratios of SFR/Ly to predict 
the fluxes of the EROs and DRGs at 1 < z < 4: we 
assume they all have the same ratio, i.e., SFR scales lin- 
early with Ly- We assume that the ratio is that of non- 
DRGs for the other galaxies at the same redshift inter- 
val. Using this very simple assumption, we find that the 
EROs+DRGs contribute 12.3 Jy deg~ 2 . This is about 
50% the background bet ween 0.5 </ < 5mJy deter- 
mined bv iKnudsenl l)2004[K which is 26 Jydeg -2 . For the 
non-EROs/DRGs at 1 < z < 4 and K < 22.5 we esti- 
mate 1.4 Jy deg -2 . 

It is interesting to speculate whether there exists 
an evolutionary link between the massive high redshift 
galaxy populations, namely the DRGs and EROs and 
the submm galaxies (SMGs) with .fs sn > 5mJy, whic h 
are found to have comparable masses (jGreve et al.1 2005'1 . 
Possibly the SMGs repr esent a short-lived ph ase with a 
duration of < 100 Myr ([Chapman et al]|2005j) . which is 
> 10 times shorter than what we find for the DRGs. The 
surface density of bright SMGs is about 0.2arcmin -2 . 
This relative to the surface density of DRGs and EROs, 
is comparable to the ratio of the duration of the phases. 

In summary, the submm detection of the DRGs 
strongly supports that these high-z, massive galaxies are 
undergoing star formation at high rates. Furthermore, 
the results presented here, based on simple assumptions, 
imply that DRGs and EROs dominate the counts around 
1 mJy. 
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